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Abstract. Applying the method of characteristics leads to wavefunctions and 
dynamic localization conditions for electrons on the one dimensional lattice under 
perpendicular time dependent electric and magnetic fields. Such conditions proceed 
again in terms of sums of products of Bessel functions of the first kind. However, 
this time one deals with both the number of magnetic flux quanta times 7r and the 
quotients between the Bloch frequency and the ones characterizing competing fields. 
Tuning the phases of time dependent modulations leads to interesting frequency mixing 
effects providing an appreciable simplification of dynamic localization conditions one 
looks for. The understanding is that proceeding in this manner, the time dependent 
superposition mentioned above gets reduced effectively to the influence of individual ac- 
fields exhibiting mixed frequency quotients. Besides pure field limits and superpositions 
between uniform electric and time dependent magnetic fields, parity and periodicity 
effects have also been discussed. 
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1. Introduction 

The dynamic localization (DL) of electrons moving on the one dimensional (ID) 
lattice under the influence of a longitudinal time dependent (TD) electric field like 
E(t) = E f(t) has received much interest since its discovery some 20 years ago [1-8]. 
This effect concerns the periodic return of the electron to the initially occupied site 
[1]. Accordingly, the mean square displacement (MSD) should remain bounded in time. 
In most cases which are of interest in practice the modulation function is periodic with 
period T, but superpositions of several ac-fields can also be considered. The DL referred 
to above should then occur under selected conditions concerning the quotients between 
the field frequencies and the Bloch frequency u>b = E ea/h, where a stands for the 
lattice spacing. Besides applications in several areas like high field and nonlinear effects 
[2,3], trapping in two level atoms [4], persistent THz emission [5,6], the generation of 
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higher harmonics [7] or the absolute negative conductance, the DL has been finally 
observed in the linear optical absorption coefficient of quantum dot superlattices [9]. 
It has also been found that the collapse of quasienergy bands is able to reflect the 
occurrence of DL [10-13]. Recent developments such as the DL of wave packets in 
barriers with TD parameters [14] , the appearance and disappearance of resonant peaks 
in I — V characteristics [15], or the influence of higher order neighbors [16], are worthy 
of being mentioned, too. 

We have to realize that the influence of the uniform magnetic field on lattice 
electrons, such as exhibited by the celebrated Harper-equation [17-19], looks quite 
interesting in many respects. The same concerns several superpositions between electric 
and magnetic fields which have been studied during course [20-22]. However, a 
systematic study of wavefunctions and of DL effects produced by superpositions for 
which both electric and magnetic fields are TD seems of having been, to the best of 
our knowledge, overlooked. We shall then use this opportunity to discuss DL effects 
provided by the longitudinal TD electric field ~E*(t) = (E(t),0, 0) working in conjunction 
with a transversal magnetic field like ~B*(t) = (0,0, B(t)), where B(t) = B g(t). To this 
aim we can start either by incorporating the electric field into the time dependent 
Harper-Hamiltonian [18], or from the decoupled limit of two parallel chains in electric 
and magnetic fields [21]. We have to realize that such systems can be converted one 
into another with the help of gauge transformations [23]. The former alternative is 
appropriate for Hall conductance studies. We shall then choose the latter alternative 
since it complies in a more suitable manner with the DL problem. This opens the way to 
the derivation of general DL conditions, but concrete realizations such as superpositions 
between uniform electric and time dependent magnetic fields, will also be discussed. 
In this context, we found that tuning the phases of TD modulations opens the way 
to interesting frequency mixing effects, which leads in turn to quickly tractable DL 
conditions proceeding again in terms of the zero's of the zero-order Bessel function of the 
first kind. The understanding is that by virtue of the phase tuning, the TD superposition 
mentioned above is able to be reduced effectively to the influence of individual ac-fields 
characterized by mixed frequency quotients. Last but not at least we shall deal with 
parity and periodicity effects. 

2. The derivation of the wavefunction 

The Hamiltonian describing the electron on the ID lattice under TD electric and 
magnetic fields specified above is given by [1,18,21] 

H = Sq I m >< m | —eaE f(t) m \ m >< m \ + 

m m 

{V^exp(-?7/2)|m + 1 >< m\ + exp(ry/2)|m >< m + 1|} (1) 

m m 
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where to is an integer ranging from — oo to oo. One has 
7 = l(t) = 2vr J- = 2ir(3 g(t) 



(2) 



where the transversal magnetic flux and the magnetic flux quantum are given by 
$ = B a 2 g(t) and <3>o = hc/e, respectively. So, there is (3q = eB a 2 /hc, which 
stands for the magnetic commensurability parameter. The constant on-site energy is 
denoted by Eq = huj , whereas V = hU is responsible for the nearest neighbor hopping 
parameter. One deals, of course, with an orthonormalized Wannier-basis for which 
< to | to' >= S mm >. We then have to account, as usual, for the TD single particle 
amplitude via 



*ff(t) >= y^C m (f) exp(-iu t) | to > 



(3) 



which leads in turn to the TD second-order discrete Schrodinger equation 

i^-C m (t) = C/exp(-i 7 /2)C m _ 1 (0+C/exp(i 7 /2)C m+1 (t)-rruv B f(t)C m (t) 
The next step is to apply the discrete Fourier-transform 
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C m (t) = — I dkexp(imk)C k (t) 



(5) 



•(4) 



where k denotes the dimensionless wave number. Then (4) becomes 

C k (t) = 



+ 2iU cos (k + + u B f(t) 



Of 



dk 



where 



(6) 



(7) 



2uo B dt uj B dt 

Now we have to remember that (4) and (6) have been discussed before when 7 = 
by resorting to the method of characteristics [1]. What then remains is to generalize 
these latter results towards incorporating 7 (t) such as given by (2). This results in the 
solution 

t 

C k {t) = exp -2iU J dt' cos f2fc(t, t') I (8) 



as it can be easily verified by direct computation. This time one has 



n k (t, t') = cos k + ^--uj B (7](t) - 7ji(0) 



(9) 



where 



dt>f(t')= V (t)- — ( 1 (t)- 1 (0)) 



2uj f 



(10) 



and 

t 

IK*) = J Mm (11) 


which is well known in the description of electric field problems. Accordingly, one gets 
faced with the normalized wavefunction 

C m (t) = exp(-im4,)J m (2U | Z(t) |) (12) 

by virtue of the well known properties of Bessel functions [24] , where 

^(f) = ^^-argZ(f) (13) 

and 

t 

Z(t) = J dt' exp(-icu B W)) ( 14 ) 


which deserve further attention. 



3. Dynamic localization effects 

Using (12) produces the MSD 

<m 2 >=^m 2 \C m (t) \ 2 =2U 2 \ Z(t) | 2 (15) 

m 

which generalizes apparently (2.7) in [1] in terms of the substitution 77(f) — > rj{t). 
Choosing as an example usual modulations like 

/(f) = cos^f) (16) 

and 

g{t) = sin(cu 2 f) (17) 

yields the characteristic function 

t 

Z(t) = J dt' exp ^-i^ sin(cuif') + info sin(w 2 t')) (18) 


which can be rewritten equivalently as 

^) = EE ex p( i M sin (^) Jn j m(7rAj) (19) 
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by virtue of expansions characterizing generating functions of Bessel functions [24], 
where 

Q m ,n = ^ (moj 2 - rvjJi) (20) 

and qj = ou B /ujj (j = 1,2). The point is to decompose Z{t) as 

Z(t) = Q 1 t + Q 2 (t) (21) 

in which Q 2 (t) oscillates with time. Having discriminated the linear term in t then 
produces the DL condition [1,13] 

Qi = (22) 

in which case the MSD remains bounded in time. On the other hand (19) shows that 
the discrimination of the linear term one looks for proceeds in terms of selected m- and 
n- values for which fl m ,n — > 0. To this aim let us assume that the frequencies lo\ and u 2 
are commensurate. This amounts to deal with quotients like 

± = "-l = *L = P ( 23) 
m uji q 2 Q 

in which P and Q are mutually prime integers. We then have to realize that the 
DL condition characterizing specifically the present superposition of TD electric and 
magnetic fields is given in terms of sums of products of Bessel functions of the first kind 
like 



CO 



Qi = Qi (<?i, ttA,) = J (<?i) Jo (tt/? ) + J^PiJpi (qi) Jqi (ttA>) = (24) 



in which 



Pi 



1 + (-l)'( p+( 3) (25) 



and P = Quj 2 juji. Using, however, g(t) = cos(u> 2 t) instead of (17) produces the DL 
condition 



Qi (?i,ttA)) 



= (26) 



Jo (qi) Jo (tt/3 ) + ^PiJpi (qi) Jqi (^Po) 
1=1 

which proceeds up to a phase factor like exp(— wr/3o), where 

Vi = exp (i^Ql) + exp (-fig*) {-lf P+Q) . (27) 

Such results indicate that contributions provided by electric and magnetic fields can be 
placed on the same footing. Of course, starting from f(t) — g(t) = cos(wit), one finds 
that the DL condition is still given by (26), but this time via P — Q — 1. 

One remarks that both (24) and (26) are sensitive to the parity of P + Q. So (24) 
becomes 

00 

Qi = QV (<Zi, ^0; 6) = Jo (qi) Jo (rc(3 )+2 ^ ^J 2Pl J 2Ql (tt/? ) = 0(28) 

1=1 
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when P + Q is an odd integer, but 

oo 

Qi = Q[ +) ^o; 6) = •/<> (gi) Jo (ttA>)+2 ^ & J P , (<Zi) Jqi W = 0(29) 

when P + Q is even. So far we have to insert = 1, but further generalizations are 
in order. Note that the presence of £j in (28) and (29) has to be understood just as an 
insertion prescription working under the sum. 

Proceeding in a similar manner one finds that (26) leads to 

Q[ +) vrA)) = Q[ +) (<Zi, ttA,; cos(ttQ//2)) = (30) 

and 

QV (<?i, *Po) = QV (<Zi, vr/5 ; cos(TrQZ)) + ^ = (31) 
when P + Q is even and odd, respectively, where 

oo 

T 1 = ry^TrA,) = 2^ sin (|(2/ + 1)Q) J (2J+1)p J (2i+1)Q (tt/3 ) .(32) 
z=i 

In the latter case the DL occurs when both real and imaginary parts of Q[ ^ (?i,7r/3 ) 
are zero, so that we have to account for twice harder computations. Such conditions are 
rather complex, so that when dealing with applications we have to resort to numerical 
studies. We can then say that looking for reasonable but tractable versions represents 
a quite desirable task. 

4. Phase tuning and frequency mixing effects 

A further point of interest is to account for the influence of phases on the DL. This 
opens the way to the derivation of reasonable DL conditions one looks for. For this 
purpose let us start from some appropriate modulations like 

/(t) = cos(wit + 0i) (33) 

and 

g{t) = cos(wit + 6 2 ) (34) 

for which u)\ = u 2 and 9j G [0,27r]. Repeating the same steps as before yields the DL 
condition 

oo 

Ql{R U Ba) = iJm(Rl)Jm(R2) = (35) 

m=— oo 

where 



i?i = qi cos(0i) + 7r/3 sin(# 2 ) 



(36) 



and 



R2 = nfio cos(# 2 ) - qi sin(^i) 



(37) 



are responsible for mixing effects concerning electric and magnetic quotients mentioned 
above. Discriminating the n = 0-term in (35) leads to the decomposition 



where 



Q l {R u R 2 ) = MR^Jq^) + AQi 



AQ 1 = 2J2 cos(7m/2) J^RJ J n (R 2 ) 

71=1 



(38) 



(39) 



Fixing, for convenience, the ^i-phase, one sees that the #2-phase can be tuned until 
R 2 = 0, in which case 



cos(6> 2 ) = -^-sin(6>i) . 

7T/J 

Accordingly, (38) becomes, interestingly enough, 
Q 1 (R 1 ,0) = J o (R 1 ) = . 



(40) 



(41) 



This shows that the celebrated DL condition derived before gets reproduced [1], now in 
terms of the mixed frequency quotient 

1/2 



Ri = qi cos(6>i) + 7i(3 



1 - 



sm 2 (6 1 ) 



(42) 



relying effectively on the influence of a single ac-field. So one gets faced with the 
simplified DL condition Ri = z n , where z = z n (n = 1,2,3, ...) stands for the root of 
Jq(z). This condition can be rewritten equivalently as 



2* + 2;s„gi sin(^i) + qf = ir 2 ffi 



(43) 



which is able to serve as an eigenvalue equation for nf3 or q 1 , respectively. However, the 
^2-phase can also be tuned so that ii/3o = 0, in which case the DL condition becomes 



Q 1 (0,R 2 ) = J (R 2 ) = 
instead of (41), where now the frequency quotient is given effectively by 



R 2 = -<7isin(6>i) + 7r/3 



1 - 



q " cosher 



nl/2 



(44) 



(45) 



One sees immediately that R 2 = z n > by virtue of (44), so that (43) is replaced by 

4 - 2z niqi cos(0i) + q\ = n 2 f3l (46) 

where, in general, n' differs from n. We then have to account for two realizations of 
the DL condition, such as indicated by (43) and (46). One realizes that dealing with 
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Ri and R 2 amounts to introduce effectively the inffuence of individual fields selected in 
accord with (42) and (45). However, we can choose Q\ such that 

sin(0i) + cos(0i) = (47) 

which also means that Q\ G [n/2,n] or Q\ G [37r/2,27r]. The corresponding DL condition 
is then given by 

R X = R 2 = z n (48) 

where by now n' = n, which represents an appreciable simplification. One realizes that 
(48) produces quantized flux values like 

A) = A)W = \ K - 2z ^ cos (^) + <?i] 1/2 (49) 

when starting from fixed gi-quotients. Such flux values exhibit unexpected limits like 

7r/3 (n) -> z n ± gi (50) 
if 6^1 — > 7r/2 and Q\ — > 37r/2, respectively. 

5. Other concrete realizations 

The superposition between a uniform electric field for which /(i) = 1 and a TD magnetic 
one deserves a little bit more attention. This time (10) becomes 

7}(t) = t - ^ («,(*) - 0(0)) (51) 

UJb 



so that 



Z(f) = ex P ( <n m*) Sm ^ mt) J TO ( qi ) J m (ttA,) (52) 

1 i m 



where now 



= ^ ( mcJ 2 - w B ) . (53) 



One sees that Z(t) contains only oscillatory contributions in so far as uj b < uj 2 . This 
means in turn that Qi = 0, so that requirements needed for the onset of DL are fulfilled 
from the very beginning. Furthermore let us assume that 

ujb = nu>2 (54) 

if ujb ^ ct>2, where n denotes a positive integer. Then the linear term gets discriminated 
again via Q m — > 0, which leads to a rather special DL condition like 

Qi = QI (52, tt/9o) =0(q 2 - 1) J^ttA,) = (55) 

where 9(x) stands for Heaviside's function. One realizes that this time the DL proceeds 
selectively in terms of zero's characterizing higher order Bessel functions for which 
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n = 1,2,3,..., in accord with (32). Complementarity, the regime is ballistic when 
Jni^Po) 7^ 0, but the DL gets restored again when the quotient ojb/^2 is either pure 
rational or irrational. 

Next let us account for pure electric and magnetic fields via (3q — > and ujb 0, 
respectively. In the first case one recovers the well-known DL condition [1] 

Qi(qi) = Jo(qi) = . (56) 

In the second case one finds 

Qi(0,ttA>) = Jo(7t/3 ) = (57) 

which means that the DL occurs whenever the number of flux quanta times 7r is a root of 
the zero order Bessel function, too. This proceeds in terms of dimensionless flux values 
centered around (3q = f3 n = z n /ir, where Jo(z n ) = 0. The approximation 

n+1 -(3 n = l (58) 

should also be mentioned, which indicates that the DL occurs periodically with unit 
dimensionless flux period. A such periodicity, which is reminiscent to Aharonov-Bohm 
oscillations, has also been remarked in superpositions between uniform magnetic and 
dc-ac electric fields [22]. 



6. Conclusions 



The influence of TD electric and magnetic fields on the DL of electrons on the ID lattice 
has been discussed systematically in terms of MSD's remaining bounded in time. The 
main result is given by the rather general formula (24). A such formula is able to exhibit 
several concrete realizations and serves as a starting point for further developments. 
Indeed, tuning the phases results in interesting frequency mixing effects which provide 
controllable DL conditions working again in terms of the zero-order Bessel function, as 
shown by (41) and (44). For this purpose one resorts to mixed frequency quotients such 
as given by (42) and (45). Such quotients proceed effectively in terms of individual ac- 
fields. Related flux quantization rules have also been derived in accord with (49), now 
by starting from fixed values of the electric frequency quotient q\ = ub/oji. Of course, 
(41) and (44) can also be discussed by starting from fixed 7T/3o-values, which leads in 
turn to the quantization of the electric field amplitude. We emphasize that such results 
can be readily generalized towards modulations for which uo\ ^ 002- We have also to 
remark that (49) can be rewritten in terms of pertinent vectors as 

T?o( n ) + tl=^n ( 59 ) 

with the understanding that 6\ stands for the angle between it n and ~~c? l . 

Pure fields limits can be readily performed, as shown by (56) and (57). One realizes 
that the first root in (57), namely the selected dimensionless magnetic flux 

fa = fa 9* 9* 0.76 (60) 

7T 
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deserves experimental verification in a close analogy with the confirmation of the DL 
concerning the "electric" quotient uob/^i — Z\ — 2.405 [9]. The superposition between 
a uniform electric field and a TD magnetic one has also been discussed. This is a rather 
special example for which the fields are accounted for in an asymmetric manner. Now 
the DL occurs in terms of (55), but when the Bloch frequency is quantized in accord 
with (54) only. Correspondingly, the electric field should be quantized itself by virtue 
of the rule 

E = 4 n) = (61) 

ea 

which can be viewed as being similar to (49). We can then say that results presented 
in this paper provide a deeper understanding of DL effects, with a special emphasis on 
the role of phase tuning effects. Applications in the design of nanoelectronic devices can 
also be invoked. 
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